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This report covers the research work performed for the period 
starting September 1991 and ending February 1992. An investigation of the 
different physical contributions in the displacement field derived from the 
variationally asymptotical analysis is performed. The analytical approach 
along with the derived displacement field and stiffness coefficients for a 
generally anisotropic thin-walled beam is presented in detail in Ref.l. A 
copy is attached in the Appendix for convenience. 

Significance of Out-of-plane Warping 

The variationally asymptotical approach does not require an a priori 
assumed displacement field and the warping function emerges as natural 
result. It follows an iterative process. The displacement function 
corresponding to the zeroth order approximation is obtained first by keeping 
the leading order terms in the energy functional. A set of successive 
corrections is added and the associated energy functional is determined. 
Corrections generating terms of the same order in the energy functional as 
previously obtained, are kept. The process is terminated when the new 
contributions generate terms of smaller order. The displacement field 
converges to the following expression: 

Uj = Ui(x)-y(s)U 2 (x)-z(s)U 3 + G(s)<p (x) 

+ gj( s)Ui(x) + + 93( s ^3( x ^ 

v 2 = U 2 (x)^+U 3 (x)^ + <p(x)r n 

v = U 2 (x)^--U 3 (x)^--<p(x)r t 

ds ds (1) 

The axial displacement is denoted by vi while v 2 and v denote the 
displacement along the tangent and normal to the cross section mid- 
surface, respectively as shown in Fig.l. The average displacement over the 
cross section along the x,y and z Cartesian coordinate system is denoted 
by Ui(x), U 2 (x) and U 3 (x), respectively. The cross sectional rotation is 

denoted by tfx). The underlined terms in Eq.(l) represent the extension and 
bending-related warping. These new terms emerges naturally in addition 
to the classical torsional-related warping G(s) 0'. They are strongly 


influenced by the material's anisotropy and vanish for materials that are 
either orthotropic or whose properties are antisymmetric relative to middle 
surface of the cross section wall. These out-of-plane warping functions 



Fig.l Coordinate system 

The contribution of out-of-plane warping was considered recently by 
Kosmatka [3 ]. Local in-plane deformations and out-of-plane warping of the 
cross section were expressed in terms of unknown functions. These 
functions wore assumed to be proportional to tbe axial strain, bending 
curvature and twist rate within the cross section and were determined 
using a finite element modeling. In our formulation, the. out-of-plane 
warping is shown to be proportional to the axial strain, bending curvature 
and twist rate. Moreover, the functions associated with each physical 
behavior are expressed in closed-form by gj(s) for the axial strain, g 2 (s) and 
g 3 (s) for the bending curvatures and G(s) for the twist rate. 

An illustration of their effect appears in Figs. 2 and 3 where the bending 
slope in a cantilevered beam is plotted along the span. The beam is 
subjected to a unit bending load at the tip and has a rectangular cross 
section with [15] 6 (Fig.2) and [30] 6 (Fig.3) layup. Two types of predictions are 
compared to the experimental results \4, 5 ]. In the first, the torsional- 
related warping is considered only while in the second the contribution of 
bending-related warping is included. Extension-related warping is 
negligible for this construction. Neglecting bending-related warping leads 
to significant errors in predictions for this case. 

Shear Deformation Contribution 

A similar behavior to the one illustrated in Figs. 2 and 3 was found in 
the theory of Ref. 5 when the shear deformation contribution is neglected. 
This may indicate that the out-of-plane warping due to bending includes 
implicitly the shear deformation contribution. In the theory of Ref.5 the 
cross section stiffness coefficients are predicted from a finite element 
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End 

Pig. 2 Bending slope in a [1516 cantilevered beam under unit tip load 
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Pig. 3 Bending slope in a [3016 cantilevered beam under unit tip load 
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simulation. The theory is not restricted to thin-walled configurations. In 
order to assess the similarity between the shear deformation contribution 
and the out-of-plane warping, the present theory and the numerical work of 
Ref. 5 are applied to the prediction of the deflection curve in a cantilevered 
beam made of graphite/epoxy material and subjected to a transverse tip 
load of 1 lb. The beam has a [15]6 layup with a rectangular cross section. 
The geometry and mechanical properties are similar to those of Ref. 5 and 
are provided in Table I. 


Table I. Cantilever Geometry and Properties 


Ply Thickness = 0.005 in 
Width = 0.923 in. 

Depth = 0.50 in. 

En = 20.6 Msi. 

E22 = E33 = 1.42 Msi. 
Gi2 =Gl3 = 0.87 Msi. 
G 23 = 0.696 Msi 
ni2 = vi3 = 0.30 
U23 = 0.34 


Figure 4 shows a similar behavior suggesting that in the present 
theory, shear deformation is implicitly accounted through bending-related 
warping. The prediction of Ref.5 are referred to as Classical when shear 
deformation is neglected. Further evidence could be provided by estimating 
the equivalent shear deformation strain in the present theory which can be 
expressed in terms of the slope of the plane that approximates the cross 
section warping. This slope is given by 

l zz 

( 2 ) 

where A and lzz denote the cross-sectional area and second moment of 

area about the z-axis, respectively. A comparison of the shear strain yxy 
over the length of the beam with the prediction of Ref. 5. is shown in Fig. 5. 
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Vertical Displacement (inches) 


The shear strain at the fixed end is 4.5924xl0" 4 based on Eq.(2) which is 
within 2 % of 4.6857xl0' 4 calculated on the basis of Ref. 5. 



End 


Pig. 4 Deflection of a [15]6 cantilevered beam under unit tip load 
Closing Remarks 

The variationally asymptotical theory developed provides a consistent 
means for including the effects of the material's anisotropy in thin-walled 
beams. Two issues have been addressed in this progress report. The first, is 
concerned with the functional form of in-plane deformation and out-of- 
plane warping contributions to the displacement field. The second, is 
concerned with the significance of shear deformation effects. 

A rigorous proof is provided for the assumed displacement field in 
Kosmatka's work [3]. Local in-plane deformations and out-of-plane 
warping of the cross section are indeed shown to be proportional to the axial 
strain, bending curvature and twist rate within the cross section. 
Moreover, their closed form functions are determined. 
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End 

Fig. 5 Shear strain in a [15]6 cantilevered beam under unit tip load 


The significance of shear deformation in the modeling of laminated 
composites was recognized in the early work of Rehfield and was followed 
by Chopra et al. by adopting a Timoshenko-type shear deformation 
formulation. The displacement field developed in the present work is shown 
to include shear deformation through the out-of-plane warping terms. A 
closed form expression for the slope of the plane that approximates the 
cross section warping is derived and shown to be within 2 % of the shear 
strain in a cantilever beam problem. 


6 



REFERENCES 


[1] , Berdichevsky, V., Armanios, E., and Badir, A., "Theory of Anisotropic 
Thin-Walled Closed Cross-Section Beams", To appear in a special issue of 
Composites Engineering , May 1992. 

[2] Armanios, E., Badir, A., and Berdichevsky, V., "Effect of damage on 
Elastically Tailored Composite Laminates", Proceedings of the AHS 
International Technical Specialists' Meeting on Rotorcraft Basic Research, 
Georgia Institute of Technology, Atlanta, Georgia, March 25-27, 1991, pp. 
48(1)-48(11). 

[3] . Kosmatka, J. B., "Extension-bend-Twist Coupling Behavior of Thin- 
walled Advanced Composite Beams with Initial Twist," Proceedings of the 
32st AIAA/ASME/AHS/ASC Structures, Structural Dynamics and 
Materials Conference, 1991, pp. 1037-1049. 

[4] . Smith, E. C., and Chopra, I., "Formulation and Evaluation of an 
Analytical Model for Composite Box-Beams," in Proceedings of the 31st 
AIAA/ASME/AHS/ASC Structures, Structural Dynamics and Materials 
Conference, 1990, pp. 759-782 

[5] . Smith, E. C., and Chopra, I., "Formulation and Evaluation of an 
Analytical Model for Composite Box-Beams,” Journal of the American 
Helicopter Society, July 1991, pp. 23-35. 

[6J. Hodges, D. H., Atilgan, A. R., Cesnik, C. S., and Fulton, M. V., "On a 
Simplified Strain Energy Function for Geometrically Nonlinear Behavior of 
Anisotropic Beams," Presented at the Seventeenth European Rotorcraft 
Forum, September 24-26, 1991, Berlin, Germany. To appear in a special 
issue of Composites Engineering, May 1992 


7 


APPENDIX 


Paper to appear in a special issue of Composites Engineering, May 1992 


Theory of Anisotropic Thin- Walled Closed 
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School of Aerospace Engineering 
Georgia Institute of Technology 
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ABSTRACT 


A variationally and asymptotically consistent theory is developed in order to denve 
the governing equations of anisotropic thin-walled beams with closed sections. e 
theory is based on an asymptotical analysis of two-dimensional shell theory. Closed- 
form expressions for the beam stiffness coefficients, stress and displacement fields ar 
provided. The influence of material anisotropy on the displacement fie.dis identified. 
A comparison of the displacement fields obtained by other analytical developments 
is performed. The stiffness coefficients and static response are also compare wi 1 
finite element predictions, closed form solutions and test data. 


INTRODUCTION 


Elastically tailored composite designs are being used to achieve favorable de f or- 
mation behavior under a given loading environment. Coupling between deformation 
modes such as extension-twist or bending-twist is created by an appropria e se ec ion 
of fiber orientation, stacking sequence and materials. The fundamental mechanism 
producing clastic tailoring in composite beams is a result of their anisotropy Sev- 
eral theories have been developed for the analysis of thin-walled anisotropic beams. 

•Professor, Associate Professor, and Graduate Research Assistant, respectively. 
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A review is provided in Hodges (1990). A basic element m the analytical model- 
ing development is the derivation of the effective stiffness coefficients and governing 
equations which allows the three-dimensional (3 D) state of stress to be recovered 
from a one-dimensional (ID) beam formulation For isotropic or ortho^pm majen- 
als this is a classical problem, which is considered m a number of text books such 
Timoshenko and Goodier (1951), Sokolnikoff (1956) Washisu (1968) Cnm<M d 
(1978), Wempner (1981), Gjelsvik (1981), Libaa and Snnmonds (1988), and Megson 

(1990). 

For generally anisotropic materials a number of ID theories have b^n developed 

by Reissner and Tsai (1972), Mansfield r “ d Sob 7 ^ 979) J 1 ™ R i9 1 91) A discussion 
(1988), Rehfield and Atilgan (1989), and Smith and Chopra (1990,1991). 

of these works is provided in the comparison section of this paper. 

The objective of this work is to develop a consistent theory for thin-walled beams 
made of anisotropic materials. The theory is an asymptotically correct first order 
approximation. The accuracy of previously developed theories is ^eje y 
in* the resulting displacement fields. A comparison of stiffness coefficients and static 
respons^with finite element predictions, closed form solutions and test data is also 

performed. 

A detailed derivation of the theory is presented first- This is followed by a num- 
mary of governing equations. Finally a comparison of results w,th previously devel- 

oped theories is provided. 

DEVELOPMENT OF THE ANALYTICAL MODEL 


Coordinate Systems 

Consider the slender thin-walled elastic cylindrical shell shown in Fig. 1. The 
length of the shell is denoted by L, its thickness by h the radius of curva urn ofthe 
middle surface by R and the maximum cross sectional dimension by d. It is assumed 

that 


d « L h « d h« R 


( 1 ) 


The shell is loaded by external forces applied to the lateral surfaces and at the 
ends It t assumed tha/thc variation of the externa, forces and proper - 

over distances of order d in the axial direction and over distances of orde 
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circumferential direction, is small. The material is anisotropic and its properties can 
vary in the direction normal to the middle surface. 

It is convenient to consider simultaneously two coordinate systems for the descrip- 
t J of Tstate of stress in thin-walled beams. The first one is teOr— 
r v and 2 shown in Fig. 1. The axial coordinate is x while y and z are associated 
with the beam cross section. The second coordinate system, is the cum ^earsys ern 
x s and f shown in Fig. 2. The circumferential coordinate s is measured g 
tangent to the middle surface in a counter-clockwise direction whereas f is measured 
Zrthe nonnal to the middle surface. A number of relationships have a ampler 
fomf when expressed in terms of curvilinear coordinates. A relationship between 
two coordinate systems can be established as follows. 

Define the position vector f of the shell middle surface as 

T = Xlx + y(s)ly + 2(s)r* 

where i x h, i z are unit vectors associated with the cartesian coordinate system X ’J 
and z. Equations y = y(s) and z = z(s) define the closed contour T in the y, z plane. 
The normal vector to the middle surface n has two nonzero components 

n = n y {s)tjj + n z (s)T z ^ 

The position vector R of an arbitrary material point can be written in the form 

R = r + £n 

Equations (2) and (3) establish the relations between the cartesian coordinates x, y, 
z and the curvilinear coordinates x, s, f. The coordinate f lies within the limits 

_Ms) < f <M£) 

2 “ q ~ 2 

The shell thickness varies along the circumferential direction and is denoted by h{s). 

The tangent vector f, the normal vector n and the projection of the position vector 
f and n are expressed in terms of the cartesian and curvilinear coordinates 


r on 


_ df _ dy_ dz~ 
t = Ts~~fc ly+ ds lz 

_ dz _ dy _ 

n = tx, ‘ = Ts h ~di' 1 
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- dy dz 

r ' = r ' t = y Ts + z Ts 

_ dz dy 

T ” = T- n = y Ts~ z Ts 

An asymptotical analysis is used to model the slender thin-walled shell as a beam 
with effective stiffnesses. The method follows an iterative process. The displacement 
function corresponding to the zeroth-order approximation is obtained first by keeping 
the leading order terms in the energy functional. A set of successive corrections is 
added to the displacement function and the associated energy functional is deter- 
mined. Corrections generating terms of the same order as previously obtained in the 
energy functional, are kept. The process is terminated when the new contributions 
do not generate any additional terms of the same order as previously obtained. 


Shell Energy Functional 


Consider in a 3 D space the prismatic shell shown in Fig. 2. A curvilinear frame x, 
s and £ is associated with the undeformed shell configuration. Values 1, 2 and 3 de- 
noting x, s, and £, respectively are assigned to the curvilinear frame. Throughout this 
section, Latin superscripts (or subscripts) run from 1 to 3, while Greek superscripts 
(or subscripts) run from 1 to 2, unless otherwise stated. 

The energy density of a 3D elastic body is a quadratic form of the strains 

U = 

The material properties are expressed by the Hookean tensor E ljkl . Following classical 
shell formulation (Koiter (1959), and Sanders (1959)) the through-the-thickness stress 
components 3 are considerably smaller than the remaining components o ere ore 


a* = 0 


( 4 ) 


The strains can be written as 


£ q 0 - 7oS + £Pa/3 


( 5 ) 


where 7a/3 and p Q p represent the in-plane strain components and the change in the 
shell middle surface curvatures , respectively. For a cylindrical shell these are relate 
to the displacement variables by 

dv\ 

' r " = a? 
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P 12 


dv\ dv 2 

2 ^=9F + ar 

dV 2 v 

122= 

d 2 u 

Pn = !h? 

_ d 2 v 1 ,dvi _ 
dsdx 4 R ds dx 

d 2 v d ,vt_, 

P22 ~ ds 2 ~ ds^R 


( 6 ) 


where v u v 7 and v represent the displacements in the axial, tangential and normal 
directions, repectively as shown in Fig. 2. These are related to the displacement 
components in cartesian coordinates by 


Vi = Ui 


where Uj, u 2 , and u 3 
respectively. 


dy dz 

V-2 = u 2~, 1" W 3'T _ 

ds ds 


dz 


dy 


v = u 'Ts- U3 Ts 


( 7 ) 


denote the displacements along the x, y and z coordinates, 


The energy density of the 2D elastic body is obtained in terms of 7a/J and p a0 by 
the following procedure. 


The 3 D energy is first minimized with respect to e a- This is equivalent to satis- 
Tying Eq. (4). The result is 

U = min U = -D a ^ 5 e a0 £ys 

Ci3 2 

where D a ^ s represents the components of the 2D moduli. The expressions for D 
are given in terms of E°^ s in the Appendix. 

The strain e Q0 from Eq. (5) is substituted into Eq. (8). After integration of the 
result over the thickness £ one obtains the energy of the shell $ per unit middle 
surface area 

2 $ = hC° l3yS 7apl'i6 + + ^CZ^PapPyS 
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where 


QaP-ti _ _ < £)a&yt > 

h 

C?™ = ^ < D a(hS Z > 

c?* = ^ < p> a/M e 2 > 

and a function of & say ot(0, between pointed brackets is defined as an integral 
through the thickness, viz., 


ot >= J 


+/»(*)/2 


h(s )/ 2 


o(0de 


( 9 ) 


For an applied external loading P„ the displacement field u t determining the 
deformed state is the stationary point of the energy functional 

(10) 


I = J $dxds - J PiUidxds 


Asymptotical Analysis of the Shell Energy Functional 


Zeroth-Order Approximation 

Let A and E be the order of displacements and stiffness coefficients C alhS , re- 
spectively- Assume that the order of the external forces is 


P~ O 



This assumption is shown later to be consistent with the equilibrium equations 
An alternative would be to assume the order of the external force as some quantity P 
Id derive the order of the displacements as PL' /Eh from an asymptotical analysis 

of the energy functional. 

For a thin-walled slender beam whose dimensions satisfy Eq. (1) the rate of change 
of the displacements along the axial direction is much smaller than their rate of change 
along the circumferential direction. That is, for each displacement component 


dvi 

dx 


« 


dvi 

ds 
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Using Eq. (6) and assuming that d is of the same order as R, the order of magnitude 
of the in-plane strains and curvatures is 

-m ~ 0 (=) 

27. j ~ 0 (^) 

7K ~ 0 (7) 

p " ~ 0 (§d 

Pn ~ 0 (^) 

~ 0 (d5) 


Since 711 and pn are much smaller than 712 , 722 and P 12 , P 22 , respectively, their 
contribution to the elastic energy is neglected. 


By keeping the leading order terms in the strain-displacement relationships, Eq. 
(6) can be written as 


2712 = 


dv\ 




dv -2 v 

722 = 57 + R 

_ (11) 

Pn 4 R ds 
_ d 2 u d ,^2 n 
P 72 -d?~ds { R ) 


The order of magnitude 
external forces is 


of the shell energy per unit area 
^ ( E A 2 h\ 


and the work 


done by 


7 


Since P,ti, << 4>, the contribution of external forces is neglected. The energy 
functional takes the form 

21 = J L f{4hC 1217 (in) 2 + 4/iC 1222 7i2722 + /iC 2222 ^) 2 + 4/i 2 C 1 1212 7i2Pi2 

+2/l 2 C, 1222 7l2P22 + 2/i 2 C 2212 722Pi2 + fc 2 Cf “ 722 P 22 
+ ^^ 2,2 (p, 2 ) 2 + y^ 222 pi2P22 + \^Cf 2 \ P 22f}dsdx (12) 

The integrand in Eq. (12) is a positive quadratic form, therefore the^mimmum of 
the functional is reached by functions v, t>i, and 112 for which ' 712 = 722 - P 12 - P 22 
0. From Eq. (11) this corresponds to 


0 

II 

(13) 

5 * + ” =0 

ds R 

(14) 

d 2 v d ( V 2 \ _ q 
ds\R) ~ 

(15) 

(14) and (15) should be single valued, i. e. 


\ds ) l J ds 

(16) 


The integral in Eq. (16) is performed along the cross sectional mid-plane closed con- 
tour T. The length of contour T is denoted by l. The bar in Eq. (1 ) an in e 
subsequent derivation denotes averaging along the closed contour r. 

Equation (13) implies that v\ is a function of x only, i.e. 

v, = Ux(x) 


(17) 


Integrate Eq. (15) to get 

t 5 = -^(x) (18) 

ds R 

where f(x) is an arbitrary function which is shown later to represent the i cross sec- 
tional rotation about the X-axis. From Eq. (1C) and (18), one obtams the relation 
between ip(x) and v 2 . 


■fM = (f) 


8 


Substitute v from Eq. (14) into Eq. (18), to get the following second-order differential 
equation for v 2 

(19) 

To solve this equation, one has to recall the relations between the radius of curvature 
R and the components y(s) and z(s) of the position vector associated with contour T 

dPz _ \_dy_ 
ds 2 Rds 

fy _ ( 20 ) 

ds 2 Rds 

It follows from Eq. (20) that g and g are solutions of the homogeneous form of Eq. 
(19) and v 2 = y?(x)r n is its particular solution. The general solution is therefore given 

= U 2 (x)-J^ + Uz(x)— + <p{x)r n (21) 

where U 2 and U 3 are arbitrary functions of x. Substitute from Eq. (21) into Eq. (14) 
to get j j 

v = U 2 (x) - U 3 (x) - <p(x)r t ( 22 ) 

Equations (17), (21) and (22) represent the curvilinear displacement field that mini- 
mizes the zeroth order approximation of the shell energy. Using Eq. (7) the curvilinear 
displacement field is written in Cartesian coordinates as 

ux = Ux(x) 

u 2 = U 2 (x) - z<p(x) 

Uz = Uz(x) + y<p(x) 

The variables Ux(x),U 2 (x) and U 3 (x) represent the average cross-sectional transla- 
tion while <p(x) the cross-sectional rotation normally referred to in beam theory as 
the torsional rotation. This displacement field corresponds to the zeroth-order ap- 
proximation and does not include bending behavior. For a centroidal coordinate 
system f/,(x), U 2 (x), U 3 (x) and <p(x) can be expressed as 

£/i(x) = 57 . 

U 2 {x) = 57 

U 3 {x) = 57 
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(u-t) 

<p(x) = ^=4 


First-Order Approximation 


A first-order approximation can be constructed by rewriting the displacement field 
in Eqs. (17), (21) and (22) in the form 


V\ = U\{x) + ttfi(s,x) 

Vi = U 2 (x)^ + t/ 3 (*)§ + <p(x)r n + w 2 (s, x) 
v = U 2 {x )^ - ^ 3 (x)^ - <p(x)r t + w(s,x ) 


(23) 


where W\ , w 2 and w can be regarded as correction functions to be determined based 
on their contributions to the energy functional. 

Substitute Eq. (23) into Eq. (6) to obtain the strains and curvatures in terms of 
the displacement corrections 

o dwi 

7n = 7 „ + 


dx 

dwo 

2712 = 27 i 2 + -zr + 27,2 , 27, 2 = 


dx 

In = 722 + 722 


dw\ 

~d 7 


dw 2 w 
172 ~~ds + ~R 


O d 2 w 

P 11 = Pu + 


o d 2 iu 
P 12 ~ P 12 + 


3 dw 2 


dsdx AR dx 


dx 2 

+ Pn 


(24) 


Pl2 — 


1 dw\ 


AR ds 


P22 — Pn + P22 


P 22 


d 2 w 
~c W 


_ A 

ds\R) 


where l° Q p and P ° Q 0 are the strains and curvatures corresponding to the zeroth-order 
approximation. These are expressed as 

7„ = C/,'(x) 

27, , = + U^(x)~ + ip'(x)r n 


7 22 — 0 
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(25) 


= !/?(*)§ - £/?(x)| - v"(x)r, 


e 1 

Pl2 " 4H 


t, = (l) s + ^ (l) s + ' p " (l)r ” 

P 22 = 0 


- <p'(*) 


The prim- in Eq. (25) denotes differentiation with respect to x. The order of 
is (4*). Among the new terms introduced by the function w { the leading ones are 
denoted by superscript * in Eq. (24). By keeping their contribution over the other 
terms, the energy functional can be represented by 

$( 711,2712 + 2712 , 722 , 0 , pn,hi) 


where terms of order ( 355 ) o r smaller such as 

/1P11712 , hPul 22 , h 2 P n p \2 , h 2 PuP 22 

h°Pi2ln , hPnln , h 2 PnPu > h 2 Pnp22 
are neglected in comparison with the following terms 

7ll7l2 , 7 h 722 , 7x2712 , 712752 

of order (#£). Similarly, the contribution of the work done by external forces, P l w i , is 
neglected since its order is (Ehg({j) in comparison with the order of the remaining 
terms in the energy functional (Ehg). Therefore in order to determine the functions 
Wi one has to minimize the functional 


j> $(7n , 27i2 + 2712 , 722, 0 , P12, P22) ds 

If the ririd body motion is suppressed the solution is unique. The terms P 12 , P 22 are 
essentiaf to the uniqueness of the solution; however, their contribution to the enerp' 
is of order (Eh^(^)) and is consequently dropped. This aspect is discussed y 
Berdichevsky and Misiura (1991) with regard to the accuracy of classical shell theory. 
The shell energy can therefore be represented by 


J = J ^$(7n,27i2 + 2712 , 722 , 0, 0,0) dsdx (26) 

It is worth noting that the bending contribution does not appear in Eq (26). That 
is to the first order approximation the shell energy corresponds to a membrane state. 
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The first variation of the energy functional is 


61 


-Ci{ 


d (2712) 


(d wA , d$ , ( du)2 ™\\ 
V ds ) + 0722 V ds R ) i 


dsdx 


(27) 


Equation (27) can be written in terms of the shear flow JV„ and hoop stress resultant 
AT m by recalling that A'n = “id W22 = The result is 

Set the first variation of the energy to zero, to obtain the following 


dN 12 

ds 

dN 22 

ds 

N22 

R 


= 0 

= 0 
= 0 


which result in (oq\ 

Nn = constant ^ ' 

“ d fe = 0 (») 

This is similar to the classical solution of constant shear flow and vanishing hoop 
stress. By setting N-n to zero the energy density is expressed in terms of 711 an 712 

0nly 2«> 1 =min24> = /l(s)(7ii) 2 + 2B(shii7i2 + c ( s )(7iii) J ( 30 ) 

722 

The variables A(s), B(s) and C(s) represent the axial, coupling and shear stiffnesses, 
mspectively. They are defined in terms of the 2D shell moduli in the Appendix. 

Equation (30) indicates that, to the first order the energy density ’ Action is 
independent of functions w, and w. That is the in-plane warping contnbution to the 

shell energy is negligible. The function w, however can be dct ® r ™ n “ i ^ * ' ( 1 

and (30) and by enforcing the condition on te, to be single valued as follows 


N _ _d*i_ _ 1 (B(s);n + C(s)-m) = constant 
d (2712) 2 - 


(31) 
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Substitute the leading terms from Eqs. (24) and (25) into Eq. (31) to get 

\bU[(x) + \c + U’M | + *'(x)r„( S ) + < 32 > 

In deriving Eq. (32) the term B^ has been neglected in comparison with 
This is possible if |£| is less or of the same order of magnitude as C. For 
case when \B\ » C additional investigation is needed. Since the elastic ene gy 
is positive definite, B 2 < AC, and B could be greater than C only if A >> C. In 
practical laminated composite designs \B\ < C, as the shear stiffness is grea er 
the extension-shear coupling. 

Equation (32) is a first-order ordinary differential equation in toi. The value of 
the constant in the right hand side of Eq. (32) can be found from the single value 
condition of function w\: 


(?B\ = 1 1 = o 

\ ds J l J ds 

The solution of Eq. (32) is determined within an arbitrary function of x. This function 
can be specified from various conditions. Each one yields aspecific interpretation o 
the variable U x . For example if m = 0 the variable Ui = v x according to Eq. [ b 
The choice of these conditions does not affect the final form of the ID beam y 
and therefore will not be specified in this formulation. The result is the following 

simple analytical solution of Eq. (32) 

(33) 


where 


w x = -yU' 2 {x) - zU’ 3 (x) + G(s)<fJ(x) + 9 x(s)U[(x) 

G(s)= /o’[ir c(T) " rnW ]‘ iT 

»w-jT 

B(s) 


b(r) - -c(t) \dr 


b{s ) = -2 


c(s) = 


1 


A e = -r\ 


(34) 


'C(s) C(s) 2 

The area enclosed by contour V is denoted by A> in Eq. (34). 

The displacement field corresponding to the first correction is obtained by sub- 
stituting Eq. (33) into Eq. (23) and dropping w 2 and w since their contribution to 
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the shell energy is negligible compared to The result referred to as first-order 
approximation is given by 

Vl = U\ (x) - y(s)ufa) - z(s)U' 3 (x ) + G(s)v'(x) + 9 i(s)U[(x) 

V2 = t/ 2 (x)^ + Uz(x)^ + <p(x)r n 

v = u ^ x) Ts ~ u ^ x) 1ls ~ V?(x)ri 


Displacement Field 

The displacement field corresponding to the next correction is found in the same 
wav A thhld correction can also be performed. However, subsequent correct.ons yteld 
only smaller terms, as shown in Badir (1992), and the displacement field converge 

to the following expression 

= U,{x) - y(s)U^x) - z(s)U' 3 (x) + G{s)<p'(x) 

+ gi{s)U[(x) + g 2 (s)U%(x) + 93(s)U 3 (x) 


V 2 = U 2 (x)^ + U 3 (x +.<p{x)r n 
v = t/ a (x)£ - U 3 (x)^ s - <f(x)r t 


(35) 


where 


9i( z ) = - J’ |«T)y(r) - jc(r) 
93 (s) = - jf’ fh( T ) z ( T ) - 


dr 


dr 


(36) 


It is seen from expressions (34) and (36) that C(s), 9i (*). 92 (a), and 93 (a) are single- 
valued functions, that is 

cm't = Gil ) = $1(0) = si (0 = 02(0) = 02(0 = 03(°) = 93 V) = 0 


The expressions for the displacements Vi, v and the first four terms in ui axe 
anatogousTo the classical theory of extension bending and torston of beams. U 
additional terms g,(s)U[, 9s(a)^' and 9s(« i" the . ; x p ress,on of u, . n Eq .^(35) 
represent warping due to axial strain and bending. eSG "T / Th strongly 

rally in addition to the classical torsional related warping G(s)<p - T y 
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influenced by the material’s anisotropy, and vanish for materials that are either or- 
thotropic or whose properties are antisymmetric relative to the shell middle surface. 
These out-of-plane warping functions were first derived by Armanios et al. (1991) for 
laminated composites. 

The contribution of out-of-plane warping was considered recently by Kosmatka 
(1991). Local in-plane deformations and out-of-plane warping of the cross section 
were expressed in terms of unknown functions. These functions were assumed to be 
proportional to the axial strain, bending curvature and twist rate within the cross 
section and were determined using a finite element modeling. In the present formula- 
tion, the out-of-plane warping is shown to be proportional to the axial strain, bending 
curvature and torsion twist rate. The functions associated with each physical behav- 
ior are expressed in closed-form by gi (s) for the axial strain, g 2 (s) and g 3 (s) for the 
bending curvatures and G{s) for the torsion twist rate. 


Strain Field 

The strain field is obtained by substituting Eq. (35) into Eq. (6) and neglecting 
terms of smaller order in the shell energy. The result is 

7ii = U[(x) - y(s)UZ(x) - z(s)US(x) 


2712 


2 At 

Iz 


c(s)(f' + 


b(s) - -c(s) 


U[ 


b(s)y(s ) - 5 c (s) 
b(s)z (s) - |c( S ) 


V'i 

US 


(37) 


722 = 0 

It is worth noting that the vanishing of hoop stress resultant in Eq, (29) and hoop 
strain in Eq. (37) should be interpreted as negligible contribution relative to other 
parameters. The longitudinal strain 711 is a linear function of y and z. This result 
was adopted as an assumption in the work of Libove (1988). 


In deriving Eq. (37), higher order terms associated with G<p" in the energy func- 
tional have been neglected in comparison with C ctp ') as shown in Badir (1992). 

This is possible if the following inequalities are satisfied 


A 

C 



« 1 



« 1 
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Constitutive Relationships 

Substitute Eq. (37) in the energy density, Eq. (30), and integrate over s to get the 
energy of ID beam theory 

I = j L - J PiUidxds ( 38 ) 

where 

= \ [C»(l/» 2 + <W) J + CssWf + c„(£7?) 2 ] 

+Cn’JW + CM u'i + C u U[U'l 

+Cnv'V'l + Cnv'U'i + C M UW . < 39 > 

Explicit expressions for the stiffness coefficients C u (t, j = 1, 4) are given in the 
Appendix. 

The constitutive relationships can be written in terms of stress resultants and kine- 
matic variables by differentiating Eq. (39) with respect to the associated ^nematic 
variable or by relating the traction T, torsional moment M x , and bending mo 
M v and M z to the shear flow and axial stress as follows 

T= J o u d£ds = j> N n ds 


M x = 0 = / j °nr n {s)d£ds = f Nx2Tn{s)ds 
= = J a u zd£ds = - f N n z{s)ds 

r. = = - j J an yd£ds = - j N ix y{s)ds 


(40) 


UU 2 * 

The shear flow N n is derived from the energy density in Eq. (31) and the axial stress 
resultant N u is given by 

= A(s) 7 n + B(s)7i2 


Nn = 


din 


and the associated axial and shear stresses are uniform through the wall thickness. 
Substitute Eq. (37) into Eqs. (31) and (41) and use Eq. (40) to get 


T ' 


Cn 

C 12 

Cl 3 

C14 


\U'x 

M x 


C\ 2 

C 22 

C 23 

C 24 

4 

<p' 

My 

► — 

Cl 3 

Cn 

C 33 

C34 


V'l 

m 2 \ 


.Cl 4 

C 2 4 

C34 

C44 . 


U'i 


( 42 ) 
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Equilibrium Equations 

The equilibrium equations can be derived by substituting the displacement field 
in Eq. (35) into the energy functional in Eq. (10) and using the pnnaple of mm, mum 
total potential energy to get 

T 1 + j> P x ds = 0 
M' x + j£ ( P z y - P y z)ds = 0 

M; + {f Pxzds)' + f Pzds = 0 . ( 43 ) 

M" + (fP x yds)' + fPyds = 0 

where P x , P y and P x are surface tractions along the x, y and 2 directions, respectively. 
One of the member of each of the following four pairs must be prescribed at the 


beam ends : 


TorUu M x or<p, M y or U 3 , and M z or U? 


(44) 


SUMMARY OF GOVERNING EQUATIONS 


The development presented in this work encompasses five equations. The first, is 
the displacement field given in Eq. (35). Its functional form was dctemined b^e 
on an asymptotical expansion of shell energy. The associated strain field is gi 
Eq (37) and the stress resultants in Eqs. (31), (40) and (41). e our , are 
constitutive relationships in Eq. (42) with the stiffness coefficients exposed 1 as m e- 
grals of material properties and cross sectional geometry in Eq. (56) of the Appendix. 
Finally the equilibrium equations and boundary conditions are given in Eq. (43) and 

(44), respectively. 

In the present development the determination of the displacement field is essential 
in obtaining accurate expressions for the beam stiffnesses. A comparison o e crive 
displacement field with results obtained by previous investigators is presented in 

following section. 


COMPARISON OF DISPLACEMENT FIELDS 
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The pioneering work of Reissner and Tsai (1972) is based on developing an exact 
solution to the governing equilibrium, compatibility and constitutive relations ips 
of shell theory. Closed as well as open cross-sections were considered. The : derived 
constitutive relationships are similar to Eq. (42). However the authors left to the 
reader the derivation of the explicit expressions for the stiffness coefficients, 
may be the reason for their work to have been overlooked. These expressions are 
important in identifying the parameters controlling the behavior and m per ormmg 
parametric design studies. Furthermore, the explicit form of the displacement fie 
helps evaluate and understand predictions of other analytical and numerical models. 

A number of assumptions were adopted in Reissner and Tsai s development re- 
garding material properties such as neglecting the coupling between m-plane strains 
and curvatures which can be significant in anisotropic materials. It is important to 
assess the influence of these assumptions on the accuracy. This has been done m the 
present work by using an asymptotical expansion of the shell energy and proving that 
the coupling and curvatures contributions to the energy are small m comparison wi 
the in-plane contribution. 

Mansfield and Sobey (1979) and Libove (1988) obtained the beano flexibilities re- 
lating the stretching, twisting and bending deformations to the applied axial load, tor- 
sional and bending moments for a special origin and axes onentation. They adopted 
the assumptions of a negligible hoop stress resultant N ss and a membrane state in 
the thin-walled beam section. Although they did not refer to the work of Reissner 
and Tsai (1972), their stiffnesses coincide for the special case outlined in Reissner and 
Tsai (1972). This special case refers to the one where the classical assumptions ot 
neglecting shear and hoop stresses and considering the shear flow to be constant is 
adopted. However, one has to carry out the details to show this fact. 

The work of Rehfield (1985) has been used in a number of composite applications. 
Rehfield’s displacement field is of the form 

' Ul = f/i(x) - y(s) [Ufa) - 2 7iy (x)] - z(s) [Ufa) - 2 7 *,(x)] + g(s,x ) 

u 2 = U 2 {x) - z(s)<p(x) ( 45 ) 

u 3 = Ufa) + y(s)<fi(x) 

where 7xr and 7xy are the transverse shear strains. The warping function g(s,x) is 

^ VCnaS p(5 ,x) = G(sV(x) (46) 


A comparison of the displacement fields in Eq. (35) and (45) shows that the warp- 
ing function in Rehfield’s formulation comprises the torsional-related contribution 
but does not include explicit terms that express the bendmg-related warping. The 
torsional warping function G(s) in Eq. (34) is different from the function in Eq. (47). 
The two expressions coincide when c = constant that is, when the wall stiffness and 
thickness are uniform along the cross section circumference. 

The torsional warping function in Eq. (47) was modified by Atilgan (1989) and 
Rehfield and Atilgan (1989) as 

' . (48) 


where 


and 


Cl 


1 




A _ Ldial 

^11 A-r, 




An 

\iA* 

An 


A 16 — 
A 66 — 


AmA?6 

An 


An J 


(49) 


(50) 


The Aij in Eq. (50) are the in-plane stiffnesses of Classical Lamination Theory 
(Jones (1975) and Vinson and Sierakowski (1987)). They are related to the modulus 
tensor by 


A„=<£ 
A 16 =< E 


1111 > , Ai 2 =<£ n22 

11,2 > , A 26 =<£ 1222 


A<y> — < E 


.2222 


^66 = < & 


1212 


A comparison of the modified torsional warping function in Eq. (48) and G{s) in 
Eq. (34) shows that they coincide for laminates with no extension-shear coupling 
£)in2 >=< £) 1222 >= 0, in Eq. (54) of the Appendix). For the case where the 
through-the-thickness contribution is neglected in Eq. (54), this reduces to A 16 — 

A%6 = 0 . 


The warping function obtained by Smith and Chopra (1990, 1991) for composite 
box-beams is identical to the expression of Rehfield and Atilgan (1989) and Atilgan 
(1989) given in Eqs. (46) and (48). 

An assessment of all the previous warping expressions can be made by checking 
whether they reduce to the exact expression for isotropic materials (see, for example, 
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Megson (1990)) 

d(s) 

with 


where fi is the shear modulus. 



[2 A e 

17W C2 ~ r " 



dr 


1 

fih(s) 


( 51 ) 


93 in Eqs. Th^t ^ ^ COnsequentIy 9u 92 and 

to <?f«W LL \V T Th t ,S the war P in g is torsion-related and reduces 
to G( s)v ^Moreover, the shear parameter c is equal to ^ and the expressions for 

<7(5) and G(s) m Eqs. (34) and (51) coincide. 

Rehfield’s warping function in Eq. (47) coincides with FY, /*, \ t, a 
properties and the thickness are uniform n i ,, Eq ’ ( 51 ) when the material 


APPLICATIONS 


entiaily asymmetric sOT^s C^/hf h y SUffneSS (CUS) “ d drcurafa - 

and Atilgan (i9S9), Hodges ^ if (1989), Rehfield 

(1990), and Smith L Ch^ra (mo'imy * * * (1990) ' Ch “ dra * * 

CUS Configuration 

This configuration produces extension-twist coupling The axial mimi- a 
in-plane stiffnesses A, B , and C given in Eo & ?•’ 00111)11112 and 

throughout the cross section Qn a vf *u ^ 53 ^ f the A PP endlx are constant 

(CUS) was adopted by Atilgan 0989? cl ™ m ^ remiaJ1 l' unifo ™ stiffness 

0989), and T 1 ^ “ 

sides are of reversed orient at ™ l' a u b , b ’ the p1 ^ la ^ u P s on opposite 
was adopted by Chid™" ‘ th f , na ™%antisymmetrie configuration 

y ^ ldJ1Q ra ci at (1990), and Smith and Chopra (1990,1991). 

Since A, B, and C are constants, the stiffness matrix in Eq. (42), for a centroidal 
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coordinate system, reduces to 


C u C l2 0 0 ' 

C 12 C 22 0 0 

0 0 C33 0 

0 0 0 Ca 4. 

The nonzero stiffness coefficients are given by 


C n = Al 


Cu ~ BA e 


C , 2 
C22 = yAs 

C 33 = A <j> z 2 ds - yr j 

£ z 2 ds 

1 

-3 

c 's> 

II 

<5 

f y z ds 


(52) 


For such a case the out-of-plane warping due to axial strain vanishes and g 1 does 
not affect the response. 


CAS Configuration 

This configuration produces bending-wist coupling. The stiffness A is constant 
throughout the cross section. For a box beam, the coupling stiffness, B m opposite 
members is of opposite sign and hence the 

ness (CAS) was adopted by Atilgan(1989), Rehiield and At. gan(19S9), Hodges e< 
0 I.(1989), and Rehfield ef ol.(1990). For a box-beam, the ply lay-ups along ^ hori- 
zontal members are mirror images, and hence the name symmetric configuration 
adopted by Chandra ef a(.(1990), and Smith and Chopra(1990,1991). The stiffness 
C in opposite members is equal. The stiffness matrix, for a centroidal system of axes, 

reduces to q q q 

0 C 22 C 2 3 0 

0 C 2 3 C33 0 

[ 0 0 0 C44J 

The nonzero stiffness coefficients are expressed by 

B 2 


Table 1: Properties of T300/5208 Graphite/Epoxy 


E n = 21-3 Msi 
E 22 = -^33 = 1-b Msi 
G 12 = G n = 0-9 Msi 
Gi 3 = 0.7 Msi 

1/12 = vis * 0-28 

v<2z — 0.5 


C 22 = 


C*23 — 


a 


Al 


2[d + a (§)] 

El A 2 
*[<+•(&)] ‘ 
A, 


C 33 = A j z 7 ds - |a - 


[ d + a (§)] 

C«-A§fis-*'* 


A e 


6C t 

Subscripts t and v denote top and vertical members, respectively. The box width 
and height are denoted by d and a, respectively. For the CAS configuration 
reference to the Cartesian coordinate system in Fig. 1, bending about the ji-axs 
coupled with torsion while extension and bending about the z-axis are decoup e 


In order to assess the accuracy of the predictions the present theory is applied to 

the box beam studied by Hodges el of. (1989)- ^ CT0SS sccUonal “ nfi 6 uratlon 15 
shown in Fig. 3 and the material properties in Table 1. 


Flexibility Coefficients 

A comparison of the flexibility coefficients Sy with the predictions from two models 
is provided in Table 2. The flexibility coefficients S XJ are obtained by inverting t e 
4x4 matrix in Eq. (42). The NABS A (Nonhomogeneous Anisotropic Beam Section 
Analysis) is a finite element model based on an extension of the work of Giavotto 
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Table 2: Comparison of Flexibility Coefficients of NABSA, TAIL and Present 
(lb,in units) 


Flexibility 

NABSA 

PRESENT 

% Diff. 

TAIL 

% Diff. 

S u x 10 s 

0.143883 

0.14491 

+0.7 

0.14491 

+0.7 

S 22 x 10 4 

0.312145 

0.32364 

+3.6 

0.32364 

+3.6 

Sn x 10 5 

-0.417841 

-0.43010 

+2.9 

-0.43010 

+ 2.9 

S 33 x 10 4 

0.183684 

0.1886 

+2.6 

0.17294 

—5.8 

S 44 x 10 5 

0.614311 

0.63429 

+3.2 

0.50157 

-18.4 


Table 3: Geometry and Mechanical Properties of Thin-Walled Beam with [+12] 4 CUS 
square cross-section 


Length = 24.0 in. 

Width = depth = 1.17 in. 
Ply thickness = 0.0075 in. 


Eu = ^22 = ^33 — 11-65 Msi 
Gy 2 = Gi3 = 0.82, Gi 3 = 0.7 Msi 
l /\2 = ^13 = 0.05, 1^23 = 0.3 


el at (1983). In this model ell possible types of warping are accounted for The 
TAIL model is based on the theory of Rehfield (1985) while neglecting the restrained 
torsional warping. The predictions of the NABSA and TAIL models are provide y 
Hodges el al( 1989). The percentage differences appearing in Table 2 are relative to 
the NABSA predictions. The present theory is in good agreement with NABSA. 
predictions show a difference ranging from +0.7 to +3.6 percent while those based 
on Rehfield’s theory (1985) range from +3.6 to -18.4 percent. 

The present theory is applied to the prediction of the tip deformation in a can- 
tilevered beam made of Graphite/Epoxy and subjected to different types ^of loa - 
ing The beam has a CUS square cross section with [+12], lay-up. The geometry 
and mechanical properties are given in Table 3. Comparisonofresultswith the 
MSC/NASTRAN finite element analysis of Nixon (1989) is provided ml able 4 Ihe 
MSC/NASTRAN analysis is based on a 2D plate model. The predictions of the 
present theory range from +1.7 to -0.7 percent difference relative to the finite elc- 
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Table 4: MSC/NASTRAN and Present Solutions for a CUS Cantilevered Beam with 
[+12] 4 Layups Subjected to Various Tip Load Cases 


Tip Load 

Tip Deformation 

% Diff. 

Axial Force (100 lb) 

Axial Force (100 lb) 

Torsional Moment (100 lb-in) 
Transverse Force (100 lb) 

NASTRAN Present 
Axial Disp. : 0.002189 in. 0.002202 in. 
Twist : 0.3178 deg. 0.32325 deg. 

Twist : 2.959 deg. 2.998 deg. 

Deflection: 1.866 in. . 1.853 in. 

+0.6 % 
+1.7 % 
+1.32 % 
-0.7 % 


Table 5: Cantilever Geometry and Properties 


Width = 0.953 in. 

Depth = 0.53 in. ■ 

Ply thickness = 0.005 in. 


E n = 20.59 Msi, Em = £33 = 1-42 Msi 
G n = G n = 0-87 Msi, G 23 = 0.7 Msi 
V \2 — ^13 = 0.42, 1/23 = 0.5 


ment results. 

For a CUS configuration, the extension-torsional response is decoupled from bend- 
in- Since C is constant and 51 does not affect the stiffness coefficients, the flexibility 
coefficients controlling extension and twist response, S n , 5 , 2 and S 22 coincide wi 
those of Atilgan (1989). and Rehfield and Atilgan (1989). As a “ n 

ial displacement and twist angle predictions comc.de. However, the lateral deflect 
under transverse load differs. The tip lateral deflectron pr^.cted using the theory f 
Rehfield (1985), and Atilgan (1989), and Rehfield and At, '| a "^ 989 ^' V 
resulting in -7.6 percentage difference compared to the NASTRAN result. 

The test data appearing in the comparisons of Figs. 4-9, are reported by Chandra 
el al (1990), and Smith and Chopra (1990, 1991). Figures 4 and 5 show the bending 
slope variation along the beam span for antisymmetric and symmetric ^cantfie 
under a 1 lb transverse tip load. The beam geometry and ^ter^ joperties are 
given in Table 5. The analytical predictions reported by Chandra et al. (1990) and 
Smith and Chopra (1990, 1991) together with results obtained on the basis o 
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analyses of Rehfield (1985), Rehfield and Atilgan (1989) .Atilgan (1989), and the 
present work axe combined in Figs. 4 and 5. Results show that the predictions of the 
present theory are the closest to the test data when compared to the other analytical 

approaches. 

The bending slope in Figs. 4 and 5 is defined in terms of the cross section rotation 
for theories including shear deformation. For the geometry and material properties 

considered, this effect is negligible as shown in Figs. 4 and 5 where the .?“3hable 
at the fixed end predicted by theories with shear deformation, is indistinguishable 
from zero. The nonzero value shown by the test data may be due to the expenmen a 
set up used to achieve clamped end conditions. 

The spanwise twist distribution of symmetric cantilevered beam with [30] e and 
[451 6 lay-ups is plotted in Figs. 6 and 7, respectively. The beams are subjected to 
a transverse tip load of 1 lb. Their dimensions and material Properties ^ given in 
Table 5. Results show that the present theory and the works of Rehfield and g 
(1989) and Atilgan (1989) are the closest to the test data. A similar behavior is 
found for the bending slope and the twist angle at the mid-span of the symmetric 
cantilevered beams appearing in Figs. 8 and 9. The beams are subjected to a tip 

torque of 1 lb-in. 


CONCLUSION 


An anisotropic thin-walled closed section beam theory has been developed based 
on an asymptotical analysis of the shell energy functional. The displacement field 
is not assumed apriori and emerges as a result of the analysis. In addition to the 
classical out-of-plane torsional warping, two new contnbutions are identified name y, 
axial strain and bending warping. A comparison of the derived governing equations 
confirms the theory developed by Reissner and Tsai. In addition, explicit closed - for ^ 
expressions for the beam stiffness coefficients, the stress and displacement fields are 
provided. The predictions of the present theory have been validated by comparison 
with finite element simulation, other closed form analyses and test data. 
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APPENDIX 


In this appendix explicit expressions for some of the relevant variables used in the 
development as well as the stillnesses C,, (i, j = 1, 4) in Eq. (42) are provide . 

The three stiffness parameters A, B and C in Eq. (30) are expressed in terms of 
the Hookean tensor E iikl as follows 


.... (< D' 1 ” >) ; 
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The 2D Young’s moduli D a ^ 6 are given by 


D+* = - E °^T~ - 
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where 
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and H„x are components of the inverse of the 2D matrix jE^ X3 • £ "E 555S "r 

Combining Eq. (34) and (53) the variables 6 and c can be written as 
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where the pointed brackets denote integration over, the thickness as defined in Eq. 

( 9 ). 

Expressions for the stiffness coefficients C 0 («, j = 1, 4) in terms of the cross 
section geometry and materials properties are as follows 
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Figure 1: Cartesian Coordinate System 


Figure 2: Curvilinear Coordinate System 


Figure 3: Beam Cross Section 


Figure 4: Bending Slope of an Anti-Symmetric [15] 6 Cantilever Under 1 lb Transverse 
Tip Load 


Figure 5: Bending Slope of a Symmetric [30], Cantilever Under 1 lb TVansverse Tip 
Load 


Figure 6: Twist of a Symmetric (30] 6 Cantilever Under 1 lb TVansverse Tip Load 


Figure 7: Twist of a Symmetric (45] 6 Cantilever Under 1 lb TVansverse Tip Load 


Figure 8: Bending slope at mid-span under unit tip torque of Symmetric lay-up 
Cantilever beams 


Figure 9: Twist at mid-span under unit tip torque of Symmetric lay-up Cantilever 
beams 
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